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We derive an analytic expression for the height correlation function of a rough surface based on 
the inverse wave scattering method of Kirchhoff theory. The expression directly relates the height 
correlation function to diffuse scattered intensity along a linear path at fixed polar angle. We test 
the solution by measuring the angular distribution of light scattered from rough silicon surfaces, 
and comparing extracted height correlation functions to those derived from atomic force microscopy 
(AFM). The results agree closely with AFM over a wider range of roughness parameters than 
previous formulations of the inverse scattering problem, while relying less on large-angle scatter data. 
Our expression thus provides an accurate analytical equation for the height correlation function of 
a wide range of surfaces based on measurements using a simple, fast experimental procedure. 

PACS numbers: 


I. INTRODUCTION 

Rough surfaces play an essential role in many physical 
phenomena including wave scatteriM , friction M, 
adhesion electrical conductivity [aU capacitance @ 
and heat transport @,0 , and in applications ranging 
from thin films to sensors [ll|, [l^ . The most direct way 
to measure morphology of rough surfaces is by scanning 
probe microscopy (SPM) O, I3l- However, probe size 
sometimes limits scanning resolution, in which case SPM 
images represent a convolution of probe tip and intrinsic 
sample geometry [IMl- Light scatter, which is faster 
and less invasive than SPM, can also measure surface 
roughness via the inverse wave scattering method [l^, 
Hi- The spatial resolution of the scattering technique 
can easily be adjusted to a desired level by tuning the 
incident wavelength A. 

One of the most important parameters for describ¬ 
ing rough surfaces is the height correlation function 
Cor{R) = {h{x + R)h{x))/a'^, where h{x) and h{x + R) 
denote height above a mean surface height {h) = 0 at 
horizontal positions x and x + R, respectively, and <j is 
the variance of h. The average is performed over x for 
fixed R. Other roughness parameters such as correla¬ 
tion length and roughness exponent can be derived from 
Cor{R) [m. 

Chandley [2^ first showed that Cor{R) could be ob¬ 
tained from a 2D Fourier transform of far-held scattered 
light intensity. However, a time-consuming series of mea¬ 
surements of scattered intensity prohles for a wide range 
of incident angles was required. In Chandley’s method, 
the height autocovariance function of the scattered wave- 
front is considered to be the same as the height autoco¬ 
variance function of the rough surface. This assumption, 
however, is only valid when wave vectors of incident and 
scattered light are close to the surface normal. Moreover, 
Chandley’s method was only practical for X ^ a. Zhao 


et al. [23j modihed Chandley’s method based on Kirch- 
hoff’s approximation and showed that Cor{R) could be 
estimated rapidly from a single measurement of scattered 
intensity along a linear detector array, for a single, ar¬ 
bitrary, hxed incident angle 9i (see Fig. 1) and an ar¬ 
bitrary wavelength A without the restriction of A ^ ct. 
Zamani et al. 2J, [2^ developed a rigorous mathematical 
foundation within Kirchhoff theory, using a saddle point 
approximation for calculating Cor{R) of very rough sur¬ 
faces from scattered light intensity measured along a spe¬ 
cial (in general curved) path along which both polar { 62 ) 
and azimuthal {(j) 2 ) scattering angles (defined in Fig. 1) 
varied. However, Zamani’s expression for Cor(R) was 
not compared directly with that derived in Re f. 12311 , nor 
with experimental results. Chakrabarti et al. [26| devel¬ 
oped an approach for obtaining Cor(R) from the angular 
dependence of the mean differential reflection coefficient 
of a one-dimensional randomly rough dielectric surface, 
also in the Kirchhoff approximation. 

In this paper, we derive a simple analytic relation be¬ 
tween Cor{R) and diffuse scattered intensity jld iff) us¬ 
ing a different mathematical approach than Ref. 23| , al¬ 
though still in the framework of inverse scattering theory 
in the Kirchhoff approximation. Our expression relates 
Cor{R) straightforwardly to (Idiff) in a simple experi¬ 
mental geometry: fixed polar scattering angle 62 , mea¬ 
sured as a function of azimuthal scattering angle 4)2 be¬ 
tween specular (^2 = 0) and an arbitrary angle far from 
specular direction. However, the derived relation be¬ 
tween Cor{R) and (Idiff) differs from that of Ref. [j^ . 
We demonstrate the accuracy of our expression by ob¬ 
taining Cor(R) from measured scattered intensity for 
surfaces whose correlation functions were obtained inde¬ 
pendently by atomic force microscopy (AFM), and com¬ 
paring the results with Cor(R) obtained from the expres¬ 
sion in Ref. [ 2 ^ . The comparison shows that the present 
expression is less reliant on accurate, large-(/) 2 , low-level 







2 


scatter data, and thus avoids unphysical fluctuations at 
R > lOrtni that can arise with the Cor{R) expression of 


Ref. j 23 | when the 4>2 measurement range is restricted. 


Moreover, our optically extracted Cor{R) agrees signifi¬ 
cantly better with AFM measurements than the Ref. [2^ 
expression for surfaces with high roughness (o’ > A). 

II. HEIGHT CORRELATION FUNCTION 
A. Derivation 

In the following derivation, for simplicity, we restrict 
our attention to self-affine fractal rough surfaces — i.e. 
surfaces for which the height difference A = {\h{xi) — 
h{x2)\) between two points xi and X2 is related to their 
separation R = |xi — a:2| by a power-law A ^ i?“, 
where a is the roughness exponent However, the 

general method is applicable to any homogeneous rough 
surface. The correlation function for self-affine surfaces 
can be written in the form Cor{R) = exp[— 
where ^ is the correlation length {Cor{R = = e~^) 

2l| — i.e. the lateral length (R) at which the correla¬ 
tion function drops to e~^ of its maximum at i? = 0 . 
The value of a can be extracted from the structure func¬ 
tion, which for a homogeneous isotropic rough surface is 
H{R) = {[h{r) — h{r + and is related to the cor¬ 

relation function by H{R) = 2(7^ [1 — Cor{R)] [Mj. For 
self-affine fractal rough surfaces the slope of the structure 
function on a log-log scale is equal to 2 a for r < ^. 

For a monochromatic incident wave 
of wave vector kmc, where r represents position, the scat¬ 
tered wave in the Kirchhoff approximation is [ 2 | 
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^ eil=lAxo+Byo+Chixo,yo)]^^^^y^ 


( 1 ) 


where 


A = sin 6*1 — sin 6*2 cos (/)2, 

B = —sin 02 sin ^2, 

C = —(cos01-I-cos02), 
a = sin0i(l — i?o) + sin02 cos())2(l-I-R q), 
b = sin 02 sin(/) 2 (l-I-Ro), 
c = cos02(l-|-R q) — cos 0 i(l — Ro). (2) 

In Eq. (I), Ro is the reflection coefficient and the integral 
is over the mean reference plane Am of the rough surface. 

The total scattered intensity includes coherent and dif¬ 
fuse parts. If the spot size is much larger than the inci¬ 
dent wavelength, the coherent intensity Ieoh appears in 
the specular direction. The root mean square of surface 
heights can be calculated from Q 



FIG. 1 : Geometry for wave scattering from a rough surface. 

where Ieoh and /q are the scattered intensity from a rough 
and smooth surface, respectively, in the specular direc¬ 
tion. For a surface with a Gaussian height probability 
distribution function of width cr, the mean diffuse inten¬ 
sity (Idiff) is related to Cor{R) by 0 

{id^ff) = 

k2 p2 r-oo 

= exp(-5) J [exp{gCor{R)) - 1] x 

JoikRVA"^ + B‘^)RdR. ( 4 ) 

Here the overbar denotes the complex conjugate of the 
scattered field F = ^(^T-^+c) depends on inci¬ 
dent and scattered angles. 


g = fc^cr^(cos 0 i -I- COS 02 )^ = k^cr'^C'^ 


( 5 ) 


is the roughness criterion, and jQ{kR^/A^ - 1 - R^) is a 
zero-order Bessel function of the first kind, which obeys 
the orthogonality relation [2^ 

^00 ^ 

/ Jo{UR)Jo{UR')UdU =-6{R-R') ( 6 ) 

Jo R 

for an arbitrary function U, where 

U = k\/A^ + B^. ( 7 ) 

for the scattering problem. To eliminate the prefactor 
in Eq. o, we define normalized diffuse intensity Id = 

(Idiff)/ (1^) e-SAM- Multiplying 1 ^ hyUJo{UR) and 
integrating yields 

pOO pOO pCiO 

/ Id{U,g)Jo{UR)UdU = / / MUR') 

Jo Jo Jo 

X {exp[5Cor(R')] - 1} MUR)R'dR'UdU 

= exp[gCor{R)] - 1, 


( 3 ) 


( 8 ) 
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where we assumed constant g and used orthogonality re¬ 
lation ([6]) to obtain the last expression. Eq. ([ 81 ) yields 
the desired analytic relation 

1 r°° 

CorlR) = - In [ / Id(U)Jo(UR)UdU -b 1 ], ( 9 ) 

9 Jo 

between Cor{R) and scattered intensity. 


B. Discussion 


In ([ 9 ]), Cor{R) is most accurately determined when 
Id(U) is known over the entire range of U values from 
zero to infinity. However, in order to derive the analytic 
relation m, it was necessary to restrict the parameter g 
to a constant value. This in turn implies, for fixed 0 i, a 
simple experimental set up with constant 62- Since U and 
g are functions of three common parameters (k, 0i,02), 
and U also a function of 02, fixing g imposes some lim¬ 
itations on U . Fixing k (or equivalently A) additionally 
limits the range of U . y /varies in magnitude 
from a minimum of zero to a maximum of 2 . Thus U is 
restricted to the range 0 —>■ 2 k. The scattering material, 
and available light sources and detectors, can further re¬ 
strict the range of k, and thus U. For example, metallic 
scatterers are highly reflective only for k < io.pl c., where 
LOp is the plasma frequency. Given such limitations, the 
question arises how accurately Cor{R) can be determined 
in a given experimental scattering configuration. 

Here we answer this question for an experimentally 
convenient scattering geometry with fixed k, a, 9 i and 
02 . Then g is fixed as required, and U varies solely as a 
function of 02- Id{U) can then be measured for selected 
fixed 02 using a conventional rectangular charge-coupled 
device (CCD) array with its axes oriented along the 02 
and 02 directions. Along the 02 direction, U then varies 
from A: | sin 6*1—sin 02 | to fc[sin^ 0 i —2 sin0i02 cos 02™°^^-l- 
sin^02]^^^ as 02 varies from 0 to Choosing 0 i = 

02 ensures that values of U down to zero are included. 
Varying 02 up to = tt includes values of U up to 

k I sin0i-|-sin02 |- FromEq. ([ 9 ]) , measuring the scattered 
light intensity vs. 02 then yields Cor{R). 


The correlation function formula in Ref.[ 22 | and its 
modification in Ref. is derived from inverse Fourier 
transform of scattered intensity and given by, 

Cor{R) = 1 + -^ ln[A J /^(fcx, fc||)e*'=ii^dfc||], (10) 


where, A:|| = A;sin(02) and k± = fccos(02) and A = 
V/ Idik±,k\\)dk\\. Here the wave vector changed to fc|| 
parallel to the surface plays the role of U in our equation, 
and exp(Ifc||r) is substituted by Jo(UR) in our equation. 
The main difference between these two equations arises 
at large R, where Eq. ( 10 ) often generates unphysical 
oscillations in Cor{R), as discussed below. 






FIG. 2: Gaussian (Top-left) and exponential (Top-right) 
height correlation functions of surfaces with correlation length 
^ = 1000 nm, rms roughness a = 100 nm, g = 1.8 calculated 
numerically from scattering of light of wavelength A = 600 
nm, incident at 61 = 50° (triangules) using the model de¬ 
scribed in the text, compared to the original height correla¬ 
tion functions (solid red curves). (Bottom-left) Deviation of 
correlation function from original one with increasing a and 
its compensation by increasing 9\,92{Bottom — right). 


C. Numerical Test 

To test our model numerically, we examine surfaces 
with two types of correlation functions (Gaussian and 
exponential) and various roughness (^, a) and scattering 
(01, 02, A, g) parameters. Substituting these parameters 
in Eq. [ 5 ] and inserting Gaussian or exponential corre¬ 
lation functions for Cor{R), we calculate Id{U) in the 
Kirchhoff approximation. We then substitute this Id{U) 
in Eq. O Fig. [ 2 ](top-left and right) shows two calculated 
Cor(R) curves, compared with the initially assumed one 
(red curves). The blue triangles show Cor{R) calculated 
for a mildly rough surface ^ = 1000 nm, cr = 100 nm 
probed with scatter parameters g = 1 . 8 , 0 i = 02 = 50 °, 
A = 600 nTO. In this case the calculated Cor[R) agrees 
very closely with the initially assumed one: coefficient of 
determination is 0.99 for both correlation functions, 
and the fit standard deviation error std = 0.004 and 0 . 02 , 
respectively, for Gaussian and exponential functions. As 
cr increases relative to A, the calculated Cor{R) increas¬ 
ingly deviates from the initial assumption (Fig. [2 bottom 
left). Reducing g by increasing 0 i, 02 compensates for 
large fccr, and recovers good agreement (Fig. [21 bottom 
right). Generally, keeping g < 8 yields good agreement, 
as expected for the Kirchhoff approximation, which is not 
accurate for g 1. 
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III. EXPERIMENTS 
A. Procedure 

To test the model experimentally, we used the light 
scattering set up shown in Fig. [S] A lens (/ = 5 cm) 
focused a He-Ne laser beam (A = 633 nm, power ^ 1 
mW, p-polarization) at 45 ° incidence angle to spot size 
Wo ~ 20/im onto unpolished back sides of commercial sil¬ 
icon wafers that were mounted on a translation stage, to 
enable convenient probing of multiple spots on each sur¬ 
face. A microscope objective of numerical aperture 0.42 
collected scattered light without polarization discrimina¬ 
tion in a cone of 25 ° half-angle around the specular di¬ 
rection. A charge-coupled device (CCD) camera placed 
in its focal plane, centered at 62 = 45 °, (j)2 = 0 °, recorded 
the intensity profile of scattered light over angular ranges 
20 ° < 6»2 < 70 °, - 25 ° < c/)2 < 25 °. All light detectable 
above noise was scattered within this cone angle, so there 
was no need to change the position of the CCD to collect 
light over a wider range. 


B. Results and Analysis 

Fig. 4 a shows a typical measured 2 D scattered intensity 
profile obtained in this configuration with 0.3 ms expo¬ 
sure. However, exposures as long as 30 ms were used 
to record low level scatter at large (j)2- Fig. 4 b shows 
a line-out along the azimuthal direction passing through 
the profile at the specular polar angle 62 = 45 °. The 
coherent specular peak and diffuse scattered light are 
both evident in this line-out. For data used in quan¬ 
titative analysis of surface roughness, it was important 
to ensure that the detector remained unsaturated over 
the entire dynamic range of scattered light intensity. To 
this end, we recorded each scattered intensity profile sev¬ 
eral times with different neutral density filters inserted 
in the path of the incident beam, so that the intense 
coherent peak and the weak tails of the diffuse profile 
were both recorded within the CCD’s linear response 
range. We then assembled the composite profile from 
the separate recordings. The shape of the scattered in¬ 
tensity profile, and extracted roughness parameters, did 
not depend significantly on polarization of the incident 
beam. From the measured linearized total scattered in¬ 
tensity profile I{4>2), we separately fitted the coherent 
central peak Icoh{4>2) and diffuse profile Idi4>2), the lat¬ 
ter from the edge of the central coherent peak out to the 
largest angle 4>2 at which scattered light was detectable 
above noise. Since samples were mounted on a transla¬ 
tion stage, we measured scattered profiles at many differ¬ 
ent spots on the surface to determine the statistical vari¬ 
ance of key roughness parameters. We then calculated a 
from Eq. (| 31 ), and Cor{R) from Eq. © or (fTUl) . Einally 
we determined a and Cor(R) independently from AEM 



FIG. 3: Experimental setup for measuring intensity profile of 
light scattered from a rough surface. Light source is a HeNe 
laser; steering mirrors enable adjustment of incident angle 0 i ; 
ND hlter ensures detector remains unsaturated; CCD camera 
is centered on the specular reflection angle, and located in the 
focal plane of a microscope objective collecting lens. 




FIG. 4: Roughness data for Si surface with cr = 175 nm 
~ A/4, (a) 2D scattered intensity profile recorded with CCD 
camera centered in the specular direction 62 = 0i = 45° and 
/>2 = 0. (b)Line-out of the profile in panel (a) along the 
azimuthal [(j> 2 ) direction, passing through the specular polar 
angle 62 = 45°, showing coherent specular peak centered at 
(j >2 = 0, and diffuse scattered profile at larger (f> 2 . (c) AFM 
image of the silicon surface, (d) correlation function for the 
silicon sample obtained from present method (red-triangle), 
AFM (black) and the method of Ref. (blue-circle). 


images, such as the one shown in Fig. 4 c, from the same 
regions of the surface from which light was scattered. 

Table I shows extracted a, ^ and a values from AFM, 
and from light scattering using the present method and 
the method of Ref. [l^, for Si surfaces with cr ^ A/ 4 . The 
uncertainties represent standard deviations of measure¬ 
ments acquired from a large number of spots on different, 
but nominally identical. Si samples. Roughness parame¬ 
ters obtained from the two light scattering methods agree 
within the stated uncertainties. The light-scatter value 
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TABLE I: Statistical surface roughness parameters obtained 
from AFM, and from light scattering using the present 
method and the method of Ref. . 



AFM 

Scattering 
(present method) 

Scattering 

(Ref. m 

method) 

cr(pm) 

0.175 ±0.002 

0.167 ±0.002 

0.167 ±0.002 

^(pm) 

6.9 ±0.2 

7.1 ±0.2 

6.8 ±0.5 

a 

0.73 ± 0.02 

0.81 ±0.02 

0.83 ±0.03 



FIG. 5: Cor{R) for Si surface with a = 657nm ~ X, derived 
with present method (red-triangle), AFM (black), and the 
method of Ref. [ 2 ^ (blue-circle). 

of cr is ^ 5% smaller, and oia ^ 10% larger, than the cor¬ 
responding AFM values, while all ^ values agree within 
measurement uncertainty. 

Fig. 4d compares Cor{R) functions obtained from 
AFM (black curve), and from light scatter analysis us¬ 
ing Eq. dni (red) and Eq. (ITOl) (blue), for the same sur¬ 
faces used for the data in Table I. Three significant differ¬ 
ences between the two light scatter analyses now emerge. 
First, the red curve, extracted from scatter profiles with 
maximum exposure 0.3 ms, did not change perceptibly 
when high-( 5 i 2 data from profiles with 30 ms exposure 
were included. In contrast, the blue curve, extracted 
from profiles with exposure up to 30 ms, changed sub¬ 
stantially when limited to data acquired with only 0.3 
ms exposure. Thus, evidently, our Eq. m relies much 
less on accurate, large-angle, low-level scatter data than 
Eq. (nnj. Second, the red curve agrees somewhat better 
with the AFM (black) curve ( 7 -^= 0 .97, standard devia¬ 
tion std=0.03) than the blue curve (r^ =0.89, std=0.11). 
Third, large fluctuations appear in the blue curve at large 
R that are not present in the other Cor{R) curves. We 
attribute these fluctuations to truncation of large fc||, 
and thus large (/) 2 ) values when evaluating the integral 
in Eq. (fT0|) . In principle, these fluctuations could be re¬ 
duced by using a wider CCD camera, or by translating 
the CCD in the azimuthal direction and using long ex¬ 


posures to acquire data over a wider ^2 range. Taken 
together, these differences demonstrate that Eq. (0 can 
yield a more accurate Cor{R) than Eq. (flOll when only 
low-exposure, low-azimuthal-angle data is available. 

Further differences between Eqs. ([9]) and cni) emerge 
when analyzing light scatter data from very rough {a > 
A) surfaces. As an example, Fig. 5 shows Cor{R) func¬ 
tions extracted from a Si surface with a — 657nm ~ A. 
The two light scatter analyses now diverge significantly 
from each other, the curve based on Eq. Q (red) agreeing 
much more closely (r^ =0.98 and std=0.03) with Cor{R) 
obtained from AFM (black curve) than the curve based 
on Eq. (flOl) (blue). Thus Eq. ([91) appears to yield more 
accurate Cor{R) than Eq. (fT0)l for very rough surfaces. 


IV. CONCLUSION 

In summary, we introduce a method for calculating 
the height correlation function Cor{R) of a homogeneous 
rough surface in the framework of Kirchhoff theory that 
yields better agreement with AFM measurements over a 
wider range of roughness parameters than previous for¬ 
mulations. Test measurements use a simple experimen¬ 
tal geometry with a single CCD camera at a fixed polar 
scattering angle that is amenable to time-resolved mea¬ 
surements. The extracted Cor{R) is free of fluctuations 
at large R that arise with previous analysis methods of 
comparable simplicity. 
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